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The intense 0.511 MeV gamma-ray line emission from the Galactic Center observed by INTE¬ 
GRAL requires a large annihilation rate of nonrelativistic positrons. If these positrons are injected 
at even mildly relativistic energies, higher-energy gamma rays will also be produced. We calculate 
the gamma-ray spectrum due to inflight annihilation and compare to the observed diffuse Galactic 
gamma-ray data. Even in a simplified but conservative treatment, we find that the positron injection 
energies must be < 3 MeV, which strongly constrains models for Galactic positron production. 

PACS numbers: 98.70.Rz, 98.70.Sa, 98.35.-a, 95.35.+d 


The central region of the Milky Way Galaxy (hereafter 
the GC) is illuminated by the annihilation of ~ 10 50 
positrons per year, producing the 0.511 MeV gamma- 
ray line flux of (1.07 ± 0.03) x 10~ 3 photons cm -2 s -1 
that is robustly detected by the INTEGRAL satellite 
0 i ! S S i 0 The angular distribution of this 
emission is aligned with the GC and is consistent with 
a 2-dimensional Gaussian of ~ 8° FWHM. There is only 
weak evidence for a disk component, in stark contrast to 
gamma-ray maps tracing nucleosynthesis (e.g., the 1.809 
MeV line from decaying 26 Al) or cosmic ray processes 
(e.g., the 1 -30 MeV continuum), which reveal a bright 
Galactic disk with several hot regions, and a much less 
prominent central region 0 j9]. Both the flux and angu¬ 
lar distribution of the 0.511 MeV radiation are difficult 
to explain Will ilium ini 

Central to resolving the origin of the positrons is the 
question of their injection energies, which ran ge u p to 
100 MeV or even higher in recent astrophysical E3 and 
exotic (requiring new particle physics) [llj models. How¬ 
ever, after energy loss (and diffusion and delay), the 
positron annihilations produce only gamma rays at or 
below 0.511 MeV, concealing their true injection energies 
(and production sites). To circumvent this, we must con¬ 
sider the relatively rare gamma-ray emission above 0.511 
MeV produced by relativistic positrons. The Internal 
Bremsstrahlung (IB) radiation associated with positron 
production (a Q.E.D. radiative correction) would conflict 
with COMPTEL and EGRET diffuse gamma-rav obser¬ 
vations unless the injection energy is < 20 MeV |12j. 

Here we consider the gamma rays produced by the In¬ 
flight Annihilation (IA) of energetic positrons with elec¬ 
trons in the interstellar medium; as noted long ago by 
Heitler, up to ~ 20% of relativistic positrons annihilate 
in flight while undergoing ionization energy loss in mat¬ 
ter [ 13 . The IA radiation as a probe of astrophysical 
positrons has also been used in Refs, EM EM EM EM EM 
, with Ref. m constraining the Galactic positrons us¬ 
ing earlier data. The new, high-quality data from INTE¬ 
GRAL allow us to model-independently study the Galac¬ 
tic positrons, normalizing the intensity and angular dis¬ 
tribution of the IA radiation to the 0.511 MeV line emis- 
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FIG. 1: Gamma-ray spectra (cusped lines, inflight annihila¬ 
tion; smooth lines, internal bremsstrahlung) from relativistic 
positrons, normalized to the 0.511 MeV flux; the injection en¬ 
ergies are 1, 3, and 10 MeV (dotted, dot-dashed and dashed). 


sion, and arriving at a very stringent limit. 

We review the energy loss and confinement of positrons 
in the GC and calculate the survival probability for rel¬ 
ativistic injected positrons to reach nonrelativistic ener¬ 
gies. We then calculate the IA gamma-ray spectrum; 
for various injection energies, Fig. ^ shows the expected 
IA and IB spectra. After comparing to Galactic diffuse 
gamma-ray data, we conclude that the IA flux is only 
compatible if the positrons are injected with energies < 3 
MeV, which limits the fraction of positrons annihilating 
in flight to < 5.5%. The generality of our constraint on 
the positron injection energy allows strong discrimination 
between models iM u • 

Positron Survival Probability. — We assume the 
positrons are injected monoenergetically at a total en¬ 
ergy Eq, and lose energy while remaining confined to the 
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GC by magnetic fields (with net displacements < 100 
pc, much less than the ~ 1 kpc size of the emission re¬ 
gion ji|). The energy loss rate for a positron of energy E 
due to ionization in a neutral hydrogen medium of num¬ 
ber density «h is (see Ref. |[oj for the more exact form 
used in our calculations) 


dE 

7.6 x 10 -26 

( nH ^ 

dx 

/3 2 

^ 0.1 cm -3 / 


(In 7 + 6 . 6 ) 


MeV 

cm 
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where 7 = E/m e = 1 / — (3 2 is the Lorentz factor and 
/3 is the velocity. For E < 100 MeV, the in-medium 
bremsstralilung and inverse Compton energy losses can 
be neglected |20}. If the medium is fully ionized, the en¬ 
ergy loss rate (Coulomb scattering instead of ionization) 
will be about 3 times larger pcj . While it has been sug¬ 
gested that the INTEGRAL positron data alone may fa¬ 
vor either a single-phase weakly-ionized warm medium [fj 
or a multi-phase medium dominated by warm neutral and 
warm ionized phases |j4) , we assume a medium of neutral 
hydrogen since direct astronomical probes suggest pfij 
a weak ionized component (~ 10%) at the GC, which 
would only mildly affect our results. 

Dirac was the first to calculate the annihilation cross 
section of positrons on electrons at rest | 22 |; it is 
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where r e is the classical electron radius. (The electrons 
at the GC can be assumed to be at rest 0, .4 0 ; when 
the electron motion must be taken into account, see 
Refs. JJ- 0 .) The energy loss is the same for every 
positron traversing the same distance (more accurately, 
the same column density). If we define the mean num¬ 
ber of positrons with energy between (E, E + dE) as 
N(E) [l5j, then the fraction of positrons annihilating as 
they travel a distance dx and lose an energy dE is 


dN(E) 

N(E) 


?lH er(-E) dx 
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where ?ih is equal to the number density of electrons 
(bound plus free). The integrated survival probability of 
positrons as they lose energy from Eq to E is [3 


N(E) ( f E ° , dE' \ 

Pe »- e= a w) =exp v" H A * E) rnm)' 

( 4 ) 

The density dependence cancels since \dE/dx\ scales with 
riH- When the positrons have lost most of their energy, 
E ~ m e ; we use P = P E 0 ^m c for this terminal survival 
probability. The argument of the exponential is small: 
for injection energies of 10 (3, 1) MeV, P differs from 
unity by ~ 11 (5.5, 1.4) %. The dominant contributions 
to IA are at increasingly lower energies, as the cross sec¬ 
tion is rising, until the energy loss rises rapidly when the 


positrons become nonrelativistic. The rate of nonrela- 
tivistic positron annihilation in the GC, N(m e ) ~ 10 50 
year -1 , has been measured by INTEGRAL. Assuming 
equilibrium of the injection and annihilation rates, the 
small IA fraction increases the required positron injec¬ 
tion rate, so that N(Eq) = N(m e )/ P. 

Nonrelativistic positrons may either directly annihi¬ 
late with an electron, producing two 0.511 MeV gamma 
rays or, due to the low temperature of the interstellar 
medium, form a positronium bound state with an elec¬ 
tron (3 23, 13] ■ Positronium annihilates to two gamma 
rays (each 0.511 MeV) 25% of the time, and to three 
gamma rays (each less than 0.511 MeV) 75% of the time. 
In the Galaxy, the relative intensities of the three-gamma 
continuum and two-gamma line emission fix the positro¬ 
nium fraction to be / = 0.967 ± 0.022 Q. The number 
of gamma rays contributing to the 0.511 MeV line per 
annihilated nonrelativistic positron is 2(1 — /) + 2/1 /4 = 
2(1 — 3//4), so that the true annihilation rate is 3.6 
times larger than would be deduced from the 0.511 MeV 
flux alone. The IA of energetic positrons produces two 
gamma rays, and so the ratio of the total flux of IA to 
0.511 MeV gamma rays is 

$ia = 2 (1 — P) = 1 1-P , , 

$0,511 2 (1 — 3//4) P 1-3//4 P ' [J> 

Since the 0.511 MeV and IA gamma rays are both emit¬ 
ted isotropically, we normalize our results to the observed 
rate and angular distribution of the 0.511 MeV data, 
eliminating the need for detailed modeling of the sources 
and the positron propagation. 

Inflight Annihilation Spectra. — The angle- 
averaged differential cross section 0 IH 0 for IA, in 
terms of the scaled gamma-ray energy k = E 1 /m e , is 

da _ nrl (-(3 + q)/(l + 7 ) + (3 + q)/fc - 1/fc 2 
dk 7 2 /3 2 \ [1 — A:/(1 + 7)] 2 

( 6 ) 

where 7(1 — /3) < 2k — 1 < 7(1 + (3). This formula 
is weighted with the gamma-ray multiplicity of 2 and 
thus integrates to twice the total cross section given in 
Eq. ©. Since the electrons are at rest, the differential 
cross section is sharply peaked at the endpoints of the k 
range. As injected positrons simultaneously lose energy 
and are annihilated, the low-energy peak remains below 
0.511 MeV, where gamma rays are accumulated, while 
the high-energy peak moves slowly down to 0.511 MeV, 
producing a long tail. Combining Eqns. 0 ® yields the 
shape of the integrated gamma-ray spectrum produced 
by the IA of positrons as they lose energy: 

d$iA _ $ 0,511 n H r E ° Ida dE 

dk ~ 1 — 3//4 P J Eo ^ E 2dk \dE/dx[ U 

The lower limit of the integral is dictated by k. due to 
energy considerations 0 . In Fig. |3 we show the IA 
gamma ray spectra for various positron injection ener¬ 
gies. The choice of plotting Ed^/dE means that the 
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variations in height directly reveal the variations in the 
numbers of gamma rays in each logarithmic energy inter¬ 
val: the area under each IA curve is proportional to twice 
the fraction of positrons annihilated in flight, as given by 
Eq. 0. The IA flux below 0.511 MeV is small compared 
to the three-gamma continuum. The gamma rays above 
0.511 MeV are spread over a broad energy range, but 
their detectability is greatly enhanced since (a) there are 
no competing signals from nonrelativistic positron anni¬ 
hilation, and (b) the Galactic diffuse gamma-ray back¬ 
ground is steeply falling with energy. We also show the 
IB gamma-ray spectrum, assuming that the positrons are 
co-produced with equal-energy electrons, which doubles 
the IB signal from positrons alone [2] . 

INTEGRAL & COMPTEL Constraints.— The 
0.511 MeV flux from the GC was measured by INTE¬ 
GRAL to be (1.07 ± 0.03) x 10 -3 photons cm -2 s _1 Q. 
The angular distribution is a two-dimensional Gaussian 
of ~ 8° FWHM (and a ~ 3.4°); this is confirmed by the 
positronium emission below 0.511 MeV [f|. This gives a 
0.511 MeV peak flux of ~ 0.048 photons cm -2 s _1 sr _1 , 
with ~ 24% (80%) of the photons coming from a circle 
of half-angle 2.5° (6°); the corresponding average fluxes 
are ~ 0.042 (0.025) photons cm -2 s _1 sr _1 . 

The flux within a circle aligned with the GC will have 
two components. First, the calculated gamma-ray spec¬ 
trum from IA (and IB) for an assumed injection energy; 
we minimize the model dependence by normalizing this 
flux and angular distribution to the observed 0.511 MeV 
data. (The IA and IB fluxes are concentrated just in this 
circle.) Second, the diffuse gamma-ray background; this 
is slowly varying (showing no excess at the GC like that 
for the 0.511 MeV line), and we normalize it from the 
measured INTEGRAL and COMPTEL data averaged in 
the inner Galactic Plane. If this combined flux at the GC 
is incompatible with the diffuse data alone, then the as¬ 
sumed positron injection energy is too large, and is disal¬ 
lowed. For monoenergetic injections, there would also be 
a sharp step in the energy spectrum inside the GC circle. 
We assume that the rates of positron production, energy 
loss, and annihilation are nearly in equilibrium when av¬ 
eraged over the propagation timescale (e.g., < 3.5 Myr 
for < 3 MeV). If positrons are not fully confined to the 
GC by magnetic fields, then the required positron pro¬ 
duction rate and the IA and IB fluxes will be even higher. 

For the diffuse flux, we used the power law 

d® ( E V 1 ' 8 o , , , , x 

— = 0.013 —— cm -2 s _1 sr _1 MeV -1 . 8 

dE \MeV J K ’ 

This reproduces the COMPTEL measurement of 0.0096 
cm -2 s _1 sr _1 between 1-3 MeV (0.0043 between 3- 
10 MeV) in the inner Galactic Plane (330° < l < 30°, 
|6| < 5°) and also agrees reasonably well with the power 
law determination below 1 MeV (with the positron¬ 
ium three-gamma continuum subtracted) from the IN¬ 
TEGRAL data alone in the region 350° < l < 10°, 
|6| < 10° [ 3 , [§, In Fig. [3 we show the diffuse spec- 



FIG. 2: The INTEGRAL and COMPTEL diffuse gamma ray 
flux measurements are shown with a black solid line, and their 
±30% uncertainties by the shaded band. For positron injec¬ 
tion energies of 1, 3, and 10 MeV (dotted, dot-dashed and 
dashed lines), the thick lines show how this would be increased 
by the inflight annihilation gamma ray flux (thin lines also in¬ 
clude the internal bremsstrahlung flux). All results are for a 
5°-diameter region at the Galactic Center. The 0.511 MeV 
line flux is not shown. 


trum, obtained by scaling the data averaged over this 
region with the solid angle of our 5°-diameter circle at 
the GC; the generous ±30% uncertainties are shown as a 
shaded band (primarily a systematic uncertainty on the 
normalization, due to subtracting detector backgrounds). 

In Fig. |21 we show how our predictions for the IA flux 
would increase the average diffuse flux in the same GC 
circle. Injection energies > 10 MeV are clearly disallowed 
since the IA flux would more than double the diffuse flux 
in this circle, giving a huge spike. Injection energies > 3 
MeV would cause a significant enhancement in the 1-3 
MeV flux in the GC circle, relative to adjoining regions. 
(Injection energies > 5 MeV give a significant enhance¬ 
ment in both the 1-3 and 3-10 MeV fluxes.) This argu¬ 
ment can be strengthened by using the COMPTEL dif¬ 
fuse skymaps directly, instead of just the averaged data 
shown in Fig. [21 We compare to the measured flux in 
strips of 5° longitude and 10° latitude @|. The 1-3 and 
3-10 MeV skymaps both show a moderate dip and peak 
structure (lower and higher flux) in strips straddling the 
GC, suggesting that its origin is in the diffuse gamma-ray 
background, and unrelated to positrons. If positrons are 
injected at > 3 MeV, then the IA flux at 1-3 MeV in 
the dip strip would be quite significant compared to the 
conventional diffuse flux there. 

Thus, even with this simplified but conservative anal¬ 
ysis, we can conclude that the positron injection energies 




4 


must be ^ 3 MeV; higher energies are very strongly ex¬ 
cluded. Including IB would strengthen our conclusions. 
Interestingly, the INTEGRAL diffuse data do show some 
excess in the 0.5-1.0 MeV flux in the GC relative to the 
adjoining regions .7]. While not strong enough to claim 
discovery of the IA flux, it is an intriguing hint, high¬ 
lighting the importance of a more sophisticated analysis, 
with optimized binning of the data in energy and angle. 

Conclusions. — The origins of the Galactic positrons 
remain mysterious. We have shown that recent high- 
quality data from INTEGRAL on the flux and angular 
distribution of the 0.511 MeV gamma-ray line enable a 
model-independent test of the origins of the positrons. 
These data allow us to accurately calculate the gamma- 
ray spectra from IA and IB produced by the positrons 
while they are relativistic, and to compare this signal to 
measurements of the diffuse ganrma-ray flux. Our results 
directly probe the positron injection energy, requiring it 
to be 3 MeV, and apply to any positron production 
mechanism. Since this is far below the energy scales sug¬ 
gested in some recent astropliysical 0 and exotic fill 
models, this is quite constraining. 

Our results complement detailed studies of nonrela- 
tivistic positron annihilation, and the prospects for a so¬ 
lution obtained with a comprehensive approach to those, 
the diffuse backgrounds, and the IA and IB fluxes are 
excellent. Since we can already easily probe injection en¬ 
ergies as small as 3 MeV, very close to the energy scale of 
nuclear beta decays from fresh nucleosynthesis products, 
this will very likely allow a first detection of the positron 
inflight annihilation signal , and INTEGRAL data appear 
to show a hint already. The angular map of positron in¬ 
jection positions given by IA and IB emission should more 
faithfully reveal the sources, in addition to the informa¬ 
tion from the injection energy scale alone. The mismatch 
between this map and the 0.511 MeV map would provide 
new insights on Galactic magnetic fields and the condi¬ 
tions of the interstellar medium. 
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